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Critical Behaviour near the Mott Metal-Insulator Transition in a 

Two-band Hubbard Model 
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The Mott metal-insulator transition in the two-band Hubbard model in infinite dimensions is 
studied by using the linearized dynamical mean-field theory. The discontinuity in the chemical 
potential for the change from hole to electron doping is calculated analytically as a function of 
the on-site Coulomb interaction U at the d-orbital and the charge-transfer energy A between 
the d- and p-orbitals. Critical behaviour of the quasiparticle weight is also obtained analytically 
as a function of U and A. The analytic results are in good agreement with the numerical results 
of the exact diagonalization method. 
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§1. Introduction 

The Mott metal- insulator transition (MIT) driven by 
the electron correlation is a fundamental problem in the 
condensed matter physics. Recently, some significant 
progress has been achieved in understanding the MIT 
by using the dynamical mean-field theory (DMFT)tl-' 
which becomes exact in the limit of infinite spatial di- 
mensions.Q In this approach, the lattice problem is 
mapped onto an effective impurity problem where a cor- 
related impurity site is embedded in an effective uncorre- 
lated medium that has to be determined self-consistently. 
To solve the effective impurity problem, several meth- 
ods have been applied including the iteratedjDerturba- 
tion theory ,EP the non-crossing approxiniation,© the pro- 
jective self-consistent method _(PSCM)p) the quantum 
Monte Carlo (QMC) methodja^ the exact diagonaliza- 
tion (ED) methoda^ and. the numerical renormalization 
group (NRG) method.lli 

In the single-band Hubbard model on the infinite di- 
mensional Bethe lattice, the Mott MIT is found to oc- 
cur as a first-order phase transition at finite temperature 
below a critical temperature Tca> Below Tc, a coexis- 
tence of the metallic and insulating solutions is found 
for the same value of the on-site CoulDpab.^iiiteractiQn 
U in the range U^T) <U < C/e2(T).00iBBBB 
At zero temperature, coexistence is also obtained for 
Uci < U < C/c2.Ba' When U increases below f/c2, the 
quasiparticle weight in the metallic solution decreases 
and finally becomes zero in the limit U -^ Uc2- When 
U decreases above C/ci, the energy gap in the insulat- 
ing solution decreases and finally becomes zero in the 
limit U —f Uci- The ground state energy in the metal- 
lic solution is lower than that in the insulating solution 
for Uci < U < Uc2- Therefore the Mott MIT occurs at 
U — Uc2 as a continuous transition at T = 0. In this 
paper we will concentrate on the Mott MIT at T = 
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and, then, we will denote the critical value Uc2 simply 

byf/e. 

The Mott MIT is observed in various id transition- 
metal compounds, which are classified into two types: 
the Mott-Hubbard (MH) type and the charge-transfer 
(CT) typeSB In the MH type such as Ti and V com- 
pounds, the Coulomb interaction U at the d-orbital is 
smaller than the CT energy A between d- and anion p- 
orbitals. In this case, the energy gap of the insulator is 
given roughly by U and a MIT occurs at a critical value 
Uc when U is varied. In the CT type such as Co, Ni and 
Cu compounds, U is larger than A. Then the energy 
gap is roughly given by A and a MIT occurs at a critical 
value Ac when A is varied. 

In this paper, we investigate the Mott MIT with both 
the MH type and the CT type. We, therefore, need to 
use the two-band Hubbard model which is characterized 
by two parameters: the on-site Coulomb interaction U 
at the d-orbital and the CT energy A between the d- and 
p-orbitals. Several authoii'iJiajsas^udied the model using 
the DMFT approachE^tiyBEl However, numerical 
problems make it difficult to study the Mott MIT for this 
model, in contrast to the single-band Hubbard model. In 
the half-filled single-band Hubbard model on the Bethe 
lattice, the chemical potential is fixed to /x = -^ because 
of the particle-hole symmetry. On the other hand, for the 
two-band Hubbard model, the Mott MIT occurs away 
from particle-hole symmetry and, then, the chemical po- 
tential has to be determined explicitly to fix the electron 
density per unit cell to be unity. This calculation con- 
sumes a lot of CPU time. Furthermore, the Mott MIT 
point is a function of several parameters, all of which 
have to be calculated, making it a difficult numerical 
problem. 

Recently, a linearized form of the DMFT has been 
developed, where the dynamicaLmepa-field equation is 
hnearized near the Mott MIT.Eiy'Ej' The linearized 
DMFT provides a simple and attractive technique to ob- 
tain approximate but analytical results for the critical 
regime of the MIT. In the single-band Hubbard model. 
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it has been demonstrated that the critical interactions of 
the MIT predicted by this approach is in very good agree- 
ment with most accurate numerical estimates^) The 
linearized DMFT has been extended to the particle-hole 
asymmetric case. This allows for a comprehensive ana- 
lytical investigation of the critical behaviour as a func- 
tioiL_nf the on-site Coulomb interaction and the dop- 
ing. Ej' The discontinuity in the chemical potential on 
changing from hole to electron doping has been calcu- 
lated analytically and is found to be in good agreement 
with the results of numerical methods, NRG as well as 
ED. Furthermore, analytic expressions for the compress- 
ibility, the quasiparticle weight, the double occupancy 
and the local spin susceptibility near half-filling, have 
been derived as functions of the interaction and the dop- 
ing. These are difficult to calculate using the numerical 
methods mentioned above. The linearized DMFT has 
also been extended to the two-band Hubbard model. Ed-* 
The phase boundary of the MIT is obtained analyti- 
cally over the whole parameter regime including the MH 
type and the CT type. The analytical result agrees well 
with the numerical result obtained from the ED method. 
However, the critical behaviour of the two-band Hubbard 
model was not considered there. 

In the present paper, we study the critical behaviour 
near the Mott MIT in the two-band Hubbard model by 
using the linearized DMFT. The model and the formu- 
lation are detailed in §2. In §3, we analytically calculate 
the discontinuity in the chemical potential for the change 
from hole to electron doping as a function of U and A 
over the hole parameter regime including the MH type 
and the CT type. The analytical results are in good 
agreement with the numerical results calculated by us- 
ing the ED method. In §4, we also analytically obtain 
the quasiparticle weight near the MIT as a function of 
U and A. The analytic results are compared with the 
ED results for several parameters. Finally, discussions 
are given in §5. 

§2. Linearized Dynamical Mean-Field Theory 

2.1 Single-band Hubbard model 

First, we consider the single-band Hubbard model, 

H = - Y^ U^j{cl^Cja + h.c.) + U'Yn,^n,i. (1) 

<i,j>,(T i 

In the limit of infinite dimensions, the self-energy be- 
comes purely site-diagonal and the local Green's func- 
tion G{z) can be given by the impurity Green's function 
of an effective single impurity Anderson model, 

+ E ^fccLcfc. + Y. Vkiflck. + cL/.), (2) 

k,(T k.a 

where £/ is the impurity level and Sk are energies of con- 
duction electrons hybridized with the impurity by Vk- In 



the nnodcl cq. (2), the non-interacting impurity Gre e n's 
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with the hybridization function. 



£k 



(4) 



includes effects of the interaction at all the sites except 
the impurity site and is determined self-consistently so 
as to satisfy the self-consistency equation. 

For simplicity, the calculations in this paper are re- 
stricted to the Bethe lattice with the connectivity q and 
the hopping tij = A=. In the limit q — oo, the self- 
consistency equation is given by 



Goiz)- 



+ H-t^G{z), 



(5) 



where fj, is the chemical potential for the original lattice 
model. Because of the particle-hole symmetry at half- 
filling, the chemical potential and the impurity level are 
set to /i = Y s-nd £/ = —^, respectively. Then, the 
self-consistency equation (ph is simply written by 

A(z) = t^G{z). (6) 

The effective impurity problem have been solved by us- 
ing various numerical methods .Ef In the metallic phase 
with intermediate interaction, the density of states is 
characterized by a three-peak structure consisting of the 
upper and the lower Hubbard bands and a quasiparticle 
peak near the Fermi level. When the system approaches 
the MIT from the metallic side at T = 0, the quasiparti- 
cle peak is found to b|e isolated from the upper and the 
lower Hubbard bands.® As U increases the width of the 
quasiparticle peak becomes narrow and finally vanishes 
in the limit U -^ Uc- r-ir-i n 

The Hnearized DMETO^yM focuses on this critical 
regime close to the MIT, where the quasiparticle peak 
is approximated by a single pole at the Fermi level, i. e. 
G{z) ~ — , near the Fermi level with a small quasipar- 
ticle weight Z ^' a,s U ^ Uc- Correspondingly, the 
hybridization function A(z) is a single-pole function, 

A(z) = — . (7) 

This represents an approximate mapping of the model 
(1) onto a two-site Anderson model,E3' 

+ e,^cK + ^E(/>- + 4/<.), (8) 

with Eg = and ef = —fj,. The hybridization strength 
V has be determined from the self-consistency equation 
which takes the simple form 



t-'z = V" 



(9) 



Now we calculate the critical value of U for the MIT at 
half-filling. In this case, the chemical potential is fixed to 



go{z)^{z-Sf-A{z))-\ 



(3) 



^ — ^ because of the particle-hob 
quasiparticle weight Z is given byE 
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^symmetry. Then the 
V (see also Appendix 
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up to second order in V. From eqs. (9) and (10), we 
obtain the critical value of the MIT in the linearized 
DMFT 



Ur = 6t. 



(11) 



The analytic result eq.(ll) is in good agreement with the 
best numericalestimates of NRG {Uc = 5.88^)0) ajid ED 
{Uc = 5.87t)El' as wen as PSCM (C/^ = 5Mt)B The 
iterated perturbation theory as well as the Gutzwiller 
approximation gives a larger critical value Uc = 6.6tcP as 
compared to the numerical approaches mentioned above. 

When we solve the self-consistency equation (9) with 
eq.(lO) by iteration, V^ increases exponentially with it- 
eration number for U < Uc and, then, the single pole 
approximation for A(z) breaks down, resulting in the 
metallic solution. For U > Uc, V^ decreases exponen- 
tially to give the self-consistent value V^ — correspond- 
ing to the insulating solution. 

For U < Uc, the chemical potential /i is continuous 
at n = 1 as a function of n. On the other hand, for 
U > Uc, n has a discontinuity at n = 1. To calculate the 
discontinuity within the linearized DMFT, we need the 
result of the quasiparticle weight Z in the |ftpneral case 



including non-symmetric case with /i ^ 2 (^^® ^^^^ 



Appendix A ) , 



Z^F{U,ti)V^ 



where 



F([/,/.) 



2^2 ^{u-^I) 2{u-^lY 



(12) 



(13) 



When U = Uc, F{U, /i) has a minimum at /i = -^ ^-nd 
t^F(U,^-^) = l which yields the critical value of the MIT 
given in eq.(ll). When U < Uc, t'^F{U,fi) > 1 for aU 
^, and the system is metallic for all n. When U > Uc, 
PF{U, ^) < 1 for /i_ < ^ < ^+ resulting in the insulat- 
ing state at n = 1, while, t^F{U,ii) > 1 for /i < /i_ or 
fj, > ijL-\. resulting in the metallic state aX n ^ 1 . Then the 
chemical potential shows a discontinuity A^ = /i+ — /x_ 
at n = 1 for the change from electron to hole doping, 
where /i± is given by the equation 



t^F{u,^,±)^l, 



(14) 



for U > Uc- By solving eq.(14) with eq.(13), we obtain 
the discontinuity in the chemical potential 



An = U il 



1 



18u2 




(15) 



where u = -g- > I. For U > Uc close to Uc, eq.(15) 
yields A^i = -^UcJ-^ - 1, while, for L/> C/c, it yields 
Afj, ~ U. The analytic result eq.(15) agrees well with th£ 
numerical resuhs of NRG and ED as weU as PSCM.B) 
We note that the discontinuity A^ calculated within the 
metallic solution is different (smaller) comuared with the 
energy gap within the insulating solution.o' In fact, Afi 
is zero for U < Uc — Uc2, while the energy gap is finite 
for U > Uci, where Ud < Uc2 as mentioned in §1. 

Finally, we study the critical behaviour near the Mott 



MIT.Eil'c3) In this case, we need the result of the quasi- 
paxtJcle weight Z up to fourth order in V , which is given 
byEJ (see also Appendix [A| ) , 
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C/2 



(16) 



in the particle-hole symmetric case with /i = 



uE> 



From 



eqs. (9) and (16), we obtain the quasiparticle weight 
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U 

TTc 



(17) 



near Uc for U < Uc- The critical property eq.(17) 
near Uc obtained from the linearized DMFT is similar 
to that from the Gutzwiller approximation which pre- 
dictsE9 Z = (1 - ■^)2. However, the coefficients C of 
the quasiparticle weight Z — C(l — ^) near Uc are dif- 
ferent, which are C = 2 within the Gutzwiller approx- 
imation and C — 18/11 within the linearized DMFT. 
The analytic result C = 18/11 seems to be still toalarge 
compared with-numerical results of NRG (C <C|-|1)I3^ and 
ED (C « 0.3)B) as weh as PSCM (C « 0.9). tl But a 
precise value of C is not obtained for the present. 

2.2 Two-band Hubbard model 

The linearized DMET has been extended to the two- 
band Hubbard model,^^ 

H = b± Y, {dlp,a + h.c.) + uY,d\^d.^dld,^ 

^ <i,j>,a i 






(18) 



This model eq.(18) is characterized by three parameters: 
the hopping integral tpd between the d- and p-orbitals, 
the on-site Coulomb interaction U at the d-orbital and 
the charge-transfer energy A = epo — e^Q between the d- 
and p-orbitals. Henceforth we set tpd — 1, unit of energy, 
and CdQ = 0, origin of energy, then, Cpo — A. 

For the model eq.(18) on the Bethe lattice with the 
connectivity q = 00, the self-consistencjj equations for 
the local Green's functions are given byED^ 



Goiz) ^ ^ z-ed-tl^Gp{z), 
Gp{z) = z — Cp — tp^Gd{z), 



(19) 
(20) 



where Gp{z) is the local Green's function for the p- 
electron and Gd{z) is that for the d-electron; e^ = 
edo - A* = -M and ep = Cpo - ^ = A - n. 

In the linearized DMFT, the two-band Hubbard model 
eq.(18) is mapped onto the two-site Anderson model 
eq.(8) with Sc ^ and ef — ed = — /i- In the limit F — > 0, 
the local Green's functions are given by Gdiz) = ^ and 
Gp{z) = -^, near the Fermi level with small weights 
Zd ^ and Zp -^ 0. Then the self-consistency equa- 
tions (19) and (20) are reduced to a simple equation 

V^. (21) 



.2 7 



To second order in V , the quasiparticle weight for the 
d-electron is given by Zd = V'^F{U, fi), and that for the 
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p-electron is given bj£3) (see also Appendix [^ 

Zj,^A{tpd,U,A,fi)V^, (22) 
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with 



A{tpd,U,A,fi) 



tldF{U,f^) 



fJ- 



2/^ 



*pd 



2(t/-M) 



(23) 



where F is defined in eq.(13). 

From eqs.(21) and (22), we have an equation to deter- 
mine the MIT point within the linearized DMFT: 



tlaAitpd,U,A,^i) 



1. 



(24) 



As mentioned in § |2.l| , V increases exponentially with 
iteration number for t^^A > 1 and, then, the single 
pole approximation for A{z) breaks down resulting in 
the metallic solution. For t^^A < 1, V^ decreases ex- 
ponentially to obtain the self-consistent value V^ = 
corresponding to the insulating solution. 

In eq.(24), A includes the chemical potential /i which 
has to be determined explicitly to obtain the critical val- 
ues of the MIT. As shown in the next paragraph, we can 
use a condition to determine /i, based on the fact that at 
the MIT point A has a minimum value as a function of 
/i. This condition gives 



—A{tpd,U,A,fi)=Q. 



(25) 



From the coupled equations (24), (25) with eq.(23), we 
obtain an analytic expression for the phase boundary 
separating the rnetallic and insulating regimes as a func- 
tion of U and A.EZP Figure. 1 shows the phase diagrain of 
the two-band Hubbard model at half-filling n — 1l3 on 
the A — U plane, where n is the electron density per unit 
cell and given by the sum of p- and d-electron densities: 
n = Up + nd- The analytic results from the linearized 
DMFT are in good agreement witLthe available numer- 
ical results from the ED methodEiJ for all values of A 
andt/O 

When the parameters, A and U, are in the metal- 
lic regime, the chemical potential /i(n) is continuous at 
n = 1 as a function of n. On the other hand, in the in- 
sulating regime, /i(ri) has a discontinuity at n = 1. Cor- 
respondingly, there are three cases in the /i dependence 
of A as below. (1) In the metallic regime, t^^^ > 1 for 
all /z resulting in the metallic solution for all n. (2) In 
the insulating regime, t^^^ < 1 for /i_ < /i < /i_|_, while, 
tpdA > 1 ioT fi < fi- or fi > fi+. Then the system is 
a Mott insulator for /i_ < /.i < jj.^, and /i shows a dis- 
continuity from /i^ to /x+ at n = 1. (3) On the phase 
boundary of the MIT, t'i^A = 1 for /i = /i(n = 1), while, 
^pd^ > 1 for /i 7^ ij.{n ~ 1). Then A has a minimum 
at /i = fi{n — 1). Therefore the equation (25) is the 
unique condition to determine the chemical potential on 
the MIT phase boundary within the linearized DMFT. 



§3. Discontinuity in the Chemical Potential 

When the parameters, A and U, are in the insulating 
regime, the chemical potential has a discontinuity A^ — 
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Fig. 1. Phase diagram of the two-band Hubbard model at T = 
and 71 = 1. Solid line is phase boundary separating the metallic 
and insulating regimes obtained from the linearized DMFT as 
a function of A and U. Closed circles are the critical values-«f 
the MIT calculated from the exact diagonalization method.L3' 
Schematic figures of the density of states are also shown for 
the Mott-Hubbard type insulator with U < A and the charge- 
transfer type insulator with U > A. 



^j^-— ji^ at n = IlS-* as mentioned in the previous section, 
where ^± is given by the equatir- 



tion 



tldA{tpd,U,A,^JL±)^l. 



(26) 



Using eq.(23) with eq.(13), we can solve eq.(26) to get 
analytic expressions for /ij- as functions of A and U , 
which yield the discontinuity A^ = fj,^ — fj,- within the 
linearized DMFT. The results of Afi are plotted as func- 
tions of A for several values of U in Fig. 2(a) and as 
functions of U for several values of A in Fig. 2(b). 

In the charge-transfer regime with U > A, the MIT 
occurs at a critical value Ac{U) when A is varied for a 
fixed U as seen in Fig. 1. In the limit U ^ oo, the critical 
value Ac and the discontinuity in the chemical potential 
A/j, are obtained by (see Appendix [Q) 



Ac = 2.08t 



pdi 



Afi = \/A2 - A2, (for A > Ac 



(27) 
(28) 



The discontinuity Afi, eq.(28), shows a square root de- 
pendence of (A — Ac) near Ac, while it is given roughly 
by A for A ^ Ac. The similar properties are also ob- 
served for finite values of U in the CT regime C/ > A as 
seen in Fig. 2(a). 

In the Mott-Hubbard regime with A > U, the MIT 
occurs at a critical value f7c(A) when U is varied for a 
fixed A as seen in Fig. 2(b). In the limit A ^ oo, the 
critical value Uc and the discontinuity in the chemical 
potential Ajj. are obtained by (see Appendix fO) 



Uc = 5.84 



A 



0, 



(29) 
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Fig. 2. Discontinuity in the chemical potential Afi obtained 
from the linearized DMFT (LDMF) as a function of A for 
[/ = 3, 4, 5, 6, 7, 8, 9, 10, oo (a) and as a function of U for A = 
3,4,5,6,7,8,9,10,00 (b). In the insets, the ED results of A^i 
are also plotted for several A at f/ = 6 (a) and for several U at 
A = 6 (b). 



A^l = U. 



(30) 



The discontinuity Afi is given roughly by U for U ^ Uc 
even in the case with finite A in the MH regime as seen 
in Fig. 2(b). 

The above mentioned features of the discontinuity in 
the chemical potential for both the CT and MH regimes 
can be clearly seen in the contour map for Afi on the 
A — t/ plane shown in Fig. 3. Note that the line with 
Afi = is equivalent to the phase boundary of the MIT 
shown in Fig. 1. 

We have also calculated A^ numerically by using the 
ED method. We solved the full DMFT equation numer- 
ically and obtained the electron density n as a function 
of /I. When fi approaches /i-|- (/i-) from above (below), 
n approaches unity from above (below) and n ~ I for 




Fig. 3. Contour map for Afi on the A — U plane calculated from 
the linearized DMFT. 



/i_ < ^ < fi^. We also confirmed that, at /i = /i±, the 
groundstate changes from singlet (/i < /i_ or /x > /i+) 
to doublet (/i_ < fi < /z+). The actual calculations were 
done for finite cluster sizes Ug up to rig = 10. The size 
dependence of Afi was very small for n^ > 6. The results 
of Afi obtained from the ED method with the cluster size 
Us — 8 are also plotted in the insets in Figs. 2(a) and (b). 
As seen in the insets in Figs. 2(a) and (b), the analytic 
results from the linearized DMFT are in good agreement 
with the numerical results from the ED method. We note 
that the value of Afi from the linearized DMFT seems to 
be larger than that from the ED for the strong-coupling 
case. This will be discussed in §5. 

§4. Critical Behavior near the Mott transition 

Finally, we discuss the critical behaviour near the Mott 
MIT at half-filling in the two-band Hubbard model. In 
this case, we need the result of the quasiparticle weight 
for the d-electron up to fourth order in V in the geneml 
case including the non-symmetric case with /i ^ ■j^j' 



(see also Appendix [~A| ) 

Zd ^ F{U, fi)V^ - G{U, fi)V\ 

where F is given in eq.(13) and 

29 24 22 



(31) 



GiU,fi) = 7r-i + 



2^4 

+ 



24 



/i2(C/-;x)2 
29 



/x(c/-m)3 2iu-fiy 



(32) 



The quasiparticle weight for the p-electron Zp is also cal- 
culated up to fourth order in V (see Appendix [^) 

Zp = Aitpd, U, A, fi)V^ - B{tpd, U, A, fi)V^ , (33) 

where A is given in eq.(23) and B is given in the Ap- 
pendix B. Substituting V^ from the self-consistency 
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equation ( pl| ) into eqs.(31) and (33), we obtain 
^ tl^A{tpd,U,A,fi)-l 



t^pdBitpd,U,A,^i) 



(34) 
(35) 



in the metallic regime close to the MIT phase boundary. 
At the MIT point with the critical values Uc, Ac and 
/Xc, eqs.(34) and (35) yield Zd — Zp — from eq.(24). 
When A or [/ decreases from the MIT point, the quasi- 
particle weight for the d-electron Zd eq.(34) and that for 
the p-elcctron Zp eq.(35) increase as 



Z,y — C/\ 



1 



z, = c-^{i- 



A 

a: 

u 



{v = d ov p) 



[v — d or p) 



(36) 



(37) 



near the MIT point at half-filling, respectively, where the 
coefficients are given by 

I _ AcF{Uc,fic) dA{tpd,Uc,A,fic) 




'^ B{tpd,U,,A,,fi,) dA 

a_ UcF{Uc,iic) dA{tpd,U,A^,iic) 



(38) 



A=A^ 



■'U 



CI 



B{tpd,Uc,Ac,fJ.c) 



dU 



, (39) 



u=u^ 
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Figures. 4(a) and (b) show the coefficients C^, Cfj, 
C\ and C\j given in eqs.(38)-(40). In the CT regime with 
(7 > A, the critical values are Ac ;$ "itpd and Uc ^ 3ipd, 
while, in the MH regime with U < A, they are Ac Z 2>tpd 
and Uc ~ 3ipd (see Fig. 1). 

In the CT regime, the MIT occurs at A = Ac for a 
fixed U Z Stpd- As A decreases below Ac for a fixed 
U, the quasiparticle weight increases as given in eq.(36). 
With increasing Uc (decreasing Ac), the coefficient C^ 
decreases due to the increasing correlation effect, while 
C^ increases because of the rapid increase in the p- 
component for the quasiparticle weight Zp/Zd — C^/C^ 
near the MIT. In the hmit Uc ^ oo {Ac -^ 2.08tpd), we Cl'"' (1 

find Ci = 1.01^^ = and C^/C^ = ClJCfj = 0.432 Unctions of A, (a) and U, (b). 
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Fig. 4. The coefficients C^ (thick sohd hne), C\ 
Hne), C? (thin solid line) and CJl (thin 



10 



Jj (thick dashed 
dashed line) of 



the quasiparticle weights for the d and p-electrons Z, 



d(p) 



-t\ and 



-'dij.) 



c\ 



d{p) 



d(p) 



(1 — ^) , respectively, as 



see Appendix |_C|) . 

In the MH regime, the MIT occurs at U = Uc for a 
fixed A < 3tpd- As U decreases below Uc for a fixed 
A, the quasiparticle weight increases as given in eq.(37). 
With increasing Ac (decreasing Uc), the coefficient C^ 
increases due to the decreasing correlation effect, while 
Ci, decreases because of the rapid decrease in the p- 
component for the quasiparticle weight Zp/Zd — C^jjCfj 
near the MIT. In the limit Ac ^ oo {Uc — > 0), we find 
Ci = 1.42, Cfj = 1.42 and C^jCl = GlJCfj = (see 
Appendix |Q) . We note that, even in the limit Ac — > oo, 
the effect of the p-band is still relevant (C^ is finite), 
because, the hopping integrals between p-p and d-d or- 
bitals were not considered in the present model eq.(18) 
and, then, the electron has to transfer between d and p 
orbitals through the hopping integral tpdS3 

We have also calculated the quasiparticle weight for 



dS(z) | _^_i 
dz 

-1 n /'~\-i 



with the local self- 
numerically by using 



the d-electron Zd = (1 -p^\z=o 

energy S(z) = Goizy^ - Gd{z) 
the ED method. The coefficients C^ and Cfj obtained 
from the linearized DMFT seem to be about five times 
larger than those from the ED method for all parame- 
ter regimes as seen in Figs. 5(a) and (b). The similar 
discrepancy has been fouad in the single-band Hubbard 
model mentioned in §2.1.ciP The trends, however, within 
the both methods are consistent with each other as fol- 
lows: (1) In the CT regime, C^ decreases with increas- 
ing U (see Fig. 5(a)). (2) In the MH regime, C^ increases 
with increasing A (see Fig. 5(b)). (3) In the intermediate 
regime C/c ~ Ac ~ 3tpd, C^ is about twice larger than 
Cf, (see Figs.5(a) and (b)). (4) C^ in the CT regime 
with U — lOtpd is almost the same as Cfj in the MH 
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(a) 



C/=3 [7=10 



ED 
LDMF 





Fig. 5. Quasiparticle weight for the d-electron Z(j calculated from 
the ED method together with that from the linearized DMFT 
as functions of A/Ac for C/ = 3 (Af ^ = 2.95) and U = 10 
(Af ° = 2.08) (a), and as functions of U/Uc for A = 3 (i7f ° = 
2.65) and A = 10 (C/f ^ = 0.60) (b). 



regime with A = lOtpd (see Figs. 5(a) and (b)). Fur- 
ther improvements in both the analytical and numerical 
methods are under way to obtain a more conclusive de- 
scription for the critical behaviour near the Mott MIT, 
which will be reported in a future publication. 

§5. Summary and Discussion 

Within the linearized DMFT, we have studied the 
Mott metal-insulator transition in the two-band Hub- 
bard model characterized by the two parameters: the d- 
site Coulomb interaction U and the d-p charge-transfer 
energy A, and obtained a unified description of the MIT 
analytically over the whole parameter regime including 
the Mott-Hubbard regime, the charge-transfer regime 
and the intermediate regime as follows: 



(1) The Mott-Hubbard regime ([/ < A) : The MIT 
occurs at a critical value Uc when U is varied for a fixed 
A at half-filling. The critical value Uc monotonically 
decreases with increasing A. When U increases above Uc 
for a fixed A, the discontinuity in the chemical potential 
AjjL increases: Ajj. oc \JU — Uc near Uc and Z\// ~ U for 
U ^ Uc- When U decreases below Uc for a fixed A, 
the quasiparticle weight for the d-electron Zd and that 
for p-electron Zp increase as Zd = C^{1 — jj-) and Zp = 
Cij{l—jj-) near Uc, respectively, where the coefficient C^ 
(Cfj) monotonically increases (decreases) with increasing 
A. 

(2) The charge-transfer regime (C/ > A) : The MIT 
occurs at Ac when A is varied for a fixed U, where 
Ac monotonically decreases with increasing U. When 
A increases above Ac for a fixed U, Afx increases as 
Z\/i ex V^ — ^c near Ac and Afi ~ A for A ^ Ac. 
When A decreases below Ac for a fixed C/, Zd and Zp 
increase as Zd = C^(l — -^) and Zp = C^(l — -^) near 
Ac, respectively, where the coefficient C^ (C'a) mono- 
tonically decreases (increases) with increasing U . 

(3) The intermediate regime (C/ ^ A) : The critical be- 
haviour near the transition as well as the phase boundary 
smoothly connects the MH type and the CT type. As U 
increases (A decreases), the p-component for the quasi- 
particle weight Zp/Zd near the MIT rapidly increases, 
which shows the change in the character of the quasipar- 
ticle from d-band like in the MH regime to p-band like 
in the CT regime. 

We have estimated the reliability of the linearized 
DMFT by comparing the present analytical results with 
the available numerical results from the exact diagonal- 
ization method. The analytical result for Zi/x neap-jthe 
MIT as well as that for the critical value of the MITO is 
in very good agreement with the numerical result within 
the ED method. In the strong coupling regime, however, 
the value of Z\/i from the linearized DMFT seems to be 
larger than that from the ED method. This may be ex- 
plained by the fact that here the chemical potential is 
very close to the edge of the Hubbard bands and/or the 
p-band, and that the effect of the bandwidth, which isjie- 
glected in the linearized DMFT, becomes importantH-* 

The analytic results for the coefficients C^ and Cfj 
from the linearized DMFT seem to be considerably larger 
than the available ED results, although a precise value 
within the purely numerical methods has not been ob- 
tained yet. When U and/or A decreases form the critical 
point of the MIT, not only the weight but also the width 
of the, quasiparticle peak is found to increase near the 
MIT.eP The width effect, which is not taken into ac- 
count in the linearized DMFT, may cause this discrep- 
ancy. Further improvements in both the analytical and 
the numerical methods are now under way. 

Effects of the hopping integrals between p-p and d-d or- 
bitals, which are not considered in the present study but 
are not negligible in actual compounds, make important 
contribution to the critical behaviour near the MIT, as 
found to be significant to determine the phase boundary 
of the MIT.EH' Furthermore, the effects of the doping, 
which is not discussed here but has been discussed in 
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the single-band Hubbard niodel,E3 make apparent dif- 
ferences between the hole doping and electron doping 
cases in the strong coupling CT regime, where the chem- 
ical potential is just below the p-band for the electron 
doping while it is just above the lower- Hubbard band for 
the hole doping. Such effects and the improvements of 
the methods will be reported in a future publication. 
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Appendix A: Two-site Anderson Model 

Here we discuss the two-site Anderson model eq.(8).Ei|) 
We assume that the conduction level is between the 
atomic /-level and the upper-Hubbard level e/ < e^ < 
Cf + U, and we set e^ = and ey = —fi. Then we define 

A = Ec - £/ = M > 0, 

Y = ef + U -ec^U -^i>0. 

Henceforth we consider the case with the small hybridiza- 
tion strength \V\ <C A, llr-and calculate all the quantities 
up to fourth order in V.c3> 

The one-electron eigenenergies and the corresponding 
eigenstates are 

y2 yi 

1/2 yi 



2 



2 



\E. 



X A3' 



where a^ = 1 — ^ -I- 3^. 

Similarly, the three-electron (one-hole) eigenenergies 
and the corresponding eigenstates are 

1/2 y4 



V 



v 



E.)^aiU + -[l~-^ cA /f /r44lO) 



Y 



where a^ = 1 - ^ -H 3^ 



Y2 



The ground state for two electron is the singlet state 
with the following eigenenergy and the eigenstate 

21/2 21/2^ 



En 



+ ef- 2V' 



1 1 

X^Y 



1 



A2 r2 



where ag ^ l-21/2( 1 + 1 ) + 41/4( 3 ^.. 

When a /, t electron is removed from the ground state 
|i?o): there are two possible final states: \E+) and \EJ). 
Correspondingly, there are two possible single-hole exci- 
tations with excitation energies, 



where 



E+-Eo = -£/ + cu2V^ - cuaV^ = -e 


1, (A-l) 


E^ - Eo = C2e2V'^ - C2e4V'^ = -£2, 


(A-2) 


3 , 2 

Cl.2 = ^ + 7, 




5 4 4 4 

r*. . 1 1 1 




"" A3 ' A2r ' Ar2 ' r3' 




1 2 

C2.2 = ^ + 7, 




3 4 4 4 




^^^^ X^ X'^Y XV^ V^ ' 





And the transition probabilities are calculated as, 
\{E+\f^\Eo)f = \- ci^2l^' + c^^aV^ 



Wl, 



\{E-\h\E 


o)^ 


= C2^2V' - 


- C2w4V'^ ; 


= W2, 






ere 

Cliu2 = 
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VA2 
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+ Ar + 
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24 


22 


24 
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12 
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C2w2 = 
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r 5 
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16 

+ A3y^ 


22 
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32 

Ay3 
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24 


C2iu2 — 


A2r2 


y4 



(A-3) 
(A-4) 



Similarly, there are two possible single-particle excita- 
tions with excitation energies. 



E^-Eo = C3,2V' 



CSeiV^ = £3, 



(A-5) 



E, 



Eq = ef + U + Ci,2V^ - Ci,4V^ = ei, (A-6) 



where 



C3e2 = 


2 
" A 
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4 
A3 


4 
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And the transition probabilities are calculated as, 



\{E_\f^\Eo)\^ = C3^2V^ ~ C3^,4V^ 



Ws, 



\{E+\ME 


0)P = 2 ~ C4r,2V^ + C4ro4V^ = W4, 


ere 

1 

C3w2 = 2 


f A 4 1 \ 
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f 24 32 22 16 5 

II II 


C3wA „ 


VX4 ' X^Y X^Y^ XY^ ' Y^ 


1 
CAw2 - 2 


(2 4 3 \ 
VX2 ' XY ' y2 j ' 


1 


( Yl 24 22 24 17 


C4iu4 „ 


VX4 X3y ^2^2 XY^ Y^ 



(A-7) 
(A-8) 



From eqs.(A-l-A-8), we obtain the /-electron Green's 
function which has four poles, G^iz) — 'Yl,i=\ -3^. In the 
limit y — > 0, high-energy poles at t\ ~ ef and £4 ~ ^f+U 
have large residues Wi « W2 ~ 5, while low-energy poles 
merge together at £2 ~ £3 ~ with small total weight 

Z = W2 + W3: 



{C2w2 + C3yj2)V'^ - {C2w4 + 03^,4)"!^'' 



(A-9) 



to fourth order in V, where 
5 4 5 

2X2 + xr ^ 2y2' 

29 24 22 



G 



24 
XY^ 



29 
2y4' 



2^4 x3r A:2y2 

Appendix B: Calculation of Zp 

In the limit V ^ 0, the local d-Green's function Gd{z) 
is approximately given by the three-pole function (see 
Appendix[A|) 

Zd 



Gd{z) 



Wi 



W4 



(B-1) 



— ei z z — 64 

where Zd ^ Z given in eq.(A-9). Substituting eq.(B-l) 
into eq.(20), Gp{z) is obtained as a four-pole function. 
In the limit V —^ 0, Gp is written by Gp{z) 



z-At 



near 



the Fermi level with the small weight Zp —> and the 
small energy Ae -^ 0. The energy Ae is calculated from 
Gp^iAe) = with eqs.(20) and (B-1): 



Ae 



2 C2«,2 + C3«;2^2 , ,2 



''pd 



Er, 



pd 



C2w4 + C3^4 

Er) 



-^{C2w2 + C3w2) + -^{C2w2 + C3u.2)^ \V^ , (B-2) 



to fourth order in V , where 

Ep = ^-^i-tl^{-^ 



1 



Bi =- 



Cle2 
2/^2 



2^ 2{U-ii), 

Clw2 , C4e2 , C4w2 



B2 = l + t 



pd 



2([/-m)2 
1 



1 

V ^ 2(C/-/i)^ 



U-fi' 



The residue Zp is calculated from Zp ^ == il'-^p ^(^)|z=/ie 
with cqs.(20), (B-1) and (B-2): 

Zp = A{tpd, U, A, ^i)V^ - B{tpd, U, A, fi)V^ , (B-3) 

to fourth order in V, where the coefficients are 

,2 C2iu2 + C3u;2 



A{tpd,U,A,fi) =t 



pd 



El 



(B-4) 



D/i 7-r A \ i2 C2tt,4 + C3^4 4 2Bi{c2w2 + 03^2) 

o[tpd, U, A, /ij = tprf -2 h tprf — 



+t 



E^ 
p 

3B2{c2w2 + C3W2) 
pd ^4 - 



E^ 
P 



(B-5) 



Appendix C: Two Limiting Cases: U ^ oo and 

A^ cx) 



In the limit U -^ cx), eqs. (24) and (25) are reduced 



to 



5tt 



2 A 



2 \ 



(C4) 



(C-2) 



respectively. By solving the coupled equations (C-1) and 
(G-2), we obtain the critical values of the MIT 



^pd 



A, ==2i 



2 2*^'^ 



1.04 i„ 



2.08 i„ 



(C-3) 



(C-4) 



Substituting eqs. (C-3) and (C-4) into eqs. (38), (39) and 
(40), we obtain 



ci = 



50(V10 - 1) 



15 + 29V10 



1.01 , 



^ _ ^5iVTo-i)i5 + AVTo)tpd 



Gfj^ 



CI 



a 



58 + 3V10 

^ = ./^-i^ 0.432, 
G^ V 5 5 



U 



0.860-^, 



in the limit U -^ oo. For A > Ac, we also solve eq.(C-l) 
to obtain /i± , which yield the discontinuity in the chem- 
ical potential Afi = fj,^ — ^_ in the limit U —^ oo: 



Afi = VA2 - A2 



(C-5) 



In the limit A —^ oo, eqs. (24) and (25) are simply 
written as 

2t^^,{5U'-2U^^ + 2^i') 

-{tl^{U-2^^) + 2{U-^^)^lAy = 0, (C-6) 

9t2^C/2 - A (lOC/3 - 22C/V -t- 6C/Ai^ - 4^^) ^ q, (C-7) 

respectively. Here we assume 

Ai = /3C/, (C-8) 
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and put eq.(C-8) into eq.(C-6), then we find 



'^ 2 19 ' 



t: 



pd 



Oi,% 7=^. 

Using eq.(C-8) with (C-9) in eq.(C-7), we have 
19 



7 



0.171, 



3a/2(341 + 95V13) 
0.419 



'^ 2 19^ 



Then the critical values of the MIT are given by 



Mc 



7 A 

U -1^ 

^^"7 A 



2.45- 



t: 



pd 



5.84^, 



(C-9) 

(C-10) 
(C-11) 

(C-12) 
(C-13) 



in the limit A -^ 00. Substituting eqs.(CT2) and (CT3) 
into eqs. (38), (39) and (40), we obtain 



^ _ 162(3051853202 + 846431785V13) 



Cfj 



348680453849 + 96706558085\/13 
15(16085875 + 4461419\/l3) 



1.42, 



114(1492608 + 413975V13) 
\/2(880\/l3 + 3173)(341 + 95\/l3)^ 



114(1492608 + 413975%/13) 



1.42, 



CI 






104831 + 29075V13 t 



'-pd 



0.895- 



t: 



pd 



0, 



18(6508 + I8O5V13) A2 ■ A2 

in the limit A -^ cx). For A = cx) and U > Uc ■ 
eq.(C-6) yields /!+ = U and /i_ == resulting in 

A^i = U. (C-14) 
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